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Abstract 

This  paper  presents  a  complete  polarization  model  of  a  solid  oxide  fuel  cell  (SOFC)  that  eliminates  the  ambiguity  of  the  suitability  of 
such  model  when  used  under  different  design  and  operating  conditions.  The  Butler- Volmer  equation  is  used  in  the  model  to  describe  the 
activation  overpotential  instead  of  using  simplified  expressions  such  as  the  Tafel  equation  and  the  linear  current-potential  equation.  In  the 
concentration  overpotential,  both  ordinary  and  Knudsen  diffusions  are  considered  to  cater  for  different  porous  electrode  designs.  Sensitivity 
tests  are  then  conducted  to  show  the  effect  of  the  thickness  of  the  respective  fuel  cell  components  on  the  drop  in  cell  voltage.  Results  show 
that  the  performance  of  an  anode-supported  fuel  cell  is  superior  to  that  using  cathode  as  the  support  under  elevated  operating  pressure  in  the 
cathode  compartment.  The  former  can  achieve  an  improved  operating  range  of  current  density  under  normal  atmospheric  conditions. 
©  2001  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

Conventional  power  stations  are,  in  general,  based  on 
vapour  power  (Rankine)  cycle  or  gas  turbine  (Brayton) 
cycle,  principles  for  power  generation.  Thus,  the  maximum 
thermal  efficiency  of  such  plants  is  restricted  by  the  Carnot 
cycle  principle.  That  is  to  say,  the  actual  thermal  efficiency 
of  a  plant  can  be  improved  either  by  increasing  the  steam/gas 
temperature  at  the  inlet  of  the  steam/gas  turbine,  or  by 
reducing  the  discharge  temperature  at  the  steam/gas  turbine 
exit.  A  high  inlet  temperature  is  normally  restricted  by 
metallurgical  considerations,  while  the  low  discharge  tem¬ 
perature  at  most  can  only  reach  the  environmental  tempera¬ 
ture  (or  the  so-called  ‘dead  state’  in  the  energy  analysis).  In 
fact,  increasing  the  inlet  temperature  of  the  steam/gas 
turbine  would  cause  a  proportional  increase  in  temperature 
at  the  steam/gas  turbine  exit  and  heat  transfer  losses,  which 
yields  only  a  marginal  improvement  thermal  efficiency 
unless  proper  recovery  of  waste  heat  is  adopted. 

In  recent  years,  the  advantages  of  gas-turbine  integrated 
fuel-cell  power  plants  with  multi-fuel  capability  that  can 
meet  commercial  performance  goals  have  been  defined  [1]. 
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Predicted  results  have  claimed  that  a  thermal  efficiency  as 
high  as  74%  lower  heating  value  (LHV)  can  be  achieved 
with  an  integrated  gas  turbine-heat  recovery  steam  generator 
(HRSG)  and  solid  oxide  fuel  cell  (SOFC)  plant.  By  contrast, 
advanced  gas  turbine-HRSG  plant  can  reach  about  62% 
LHV  [2].  All  models  of  fuel  cell  in  the  simulation  code  for 
predicting  the  plant  efficiency  rely  extensive  experimental 
data  to  characterize  the  polarization  curve  of  the  fuel  cell 
sub-system,  given  that  all  cells  have  the  same  characteristics. 
It  is  well-known  that  the  polarization  curve  is  a  function  of 
both  the  operating  conditions  and  the  component  design  of 
the  fuel  cell.  The  operating  conditions  include  the  working 
temperature,  partial  pressures  of  the  fuel  and  oxidant  and 
their  utilization  rate,  and  the  water  concentration  in  the 
anode  compartment.  The  design  specifications  include  the 
porosity,  tortuosity  and  thickness  of  the  anode  and  cathode 
(responsible  for  concentration  overpotential),  the  thickness 
of  the  electrolyte  (responsible  for  ohmic  overpotential),  and 
the  electrode/electrolyte  interface  (responsible  for  activation 
overpotential).  A  complete  polarization  model,  if  available, 
will  be  extremely  useful  to  characterize  the  operation  of  the 
fuel  cell  sub-system  provided  that  the  design  specifications, 
as  mentioned  above,  are  given.  In  addition,  a  generalized 
model  would  replicate  the  fuel  cell  performance  reliably 
under  off-design  operating  conditions  such  as  when  the 
current  densities  are  away  from  those  at  the  operating  power. 
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Nomenclature 

C 

concentration  (mol  m  3) 

D 

diffusion  coefficient  (m2  s  ') 

e 

energy  of  molecular  interaction  (ergs) 

E 

voltage  (V) 

F 

Faraday  constant  (=96485  C  mol-1) 

8 

molar  Gibbs  free  energy  (J  mol-1) 

G 

Gibbs  free  energy  (J) 

h 

molar  enthalpy  (J  mol-1) 

i 

current  density  (A  m-2) 

io 

exchange  current  density  (A  m-2) 

Ga 

anode  exchange  current  density  (A  m-2) 

Gc 

cathode  exchange  current  density  (A  m-2) 

J 

flux  (mol  s-1) 

K 

equilibrium  constant 

4 

anode  thickness  (m) 

4 

cathode  thickness  (m) 

4 

electrolyte  thickness  (m) 

M 

molecular  mass 

n 

number  of  moles 

ne 

electrons  transferred  per  reaction 

P 

pressure  (Pa) 

Q 

heat  (J) 

Q 

heat  transfer  rate  (J  s-1) 

r 

average  pore  radius  (m) 

R 

universal  gas  constant  (=8.314  J  mol-1  K-1) 

R 

area-specific  resistance  (Q  m-2) 

Re 

electrolyte  area-specific  ohmic  resistance 

(Q  m-2) 

Rc  one 

area-specific  concentration  resistance  (Q  m-2) 

S 

molar  entropy  (J  mol-1  K  ') 

SA 

surface  area  of  the  porous  solid  (m2  kg-1) 

S,  a 

sensitivity  of  cell  voltage  due  to  the  change  of 

anode  thickness 

Sk 

sensitivity  of  cell  voltage  due  to  the  change  of 

cathode  thickness 

Sk 

sensitivity  of  cell  voltage  due  to  the  change  of 

electrolyte  thickness 

W 

rate  of  work  (J  s-1) 

Wei 

reversible  electrical  work  (J) 

T 

temperature  (K) 

V 

velocity  of  the  gas  molecule  (m  s-1) 

X 

mole  fraction 

z 

number  of  electrons  participating  in  the  reaction 

Greek  Letters 

P 

transfer  coefficient 

£ 

porosity 

n 

polarization  (V) 

V  Act 

activation  polarization  (V) 

Pc  one 

concentration  polarization  (V) 

>7  Ohm 

ohmic  polarization  (V) 

P 

resistivity  (Q  m) 

Pb 

bulk  density  of  the  solid  particle  (kg  m-3) 

er 

collision  diameter  (A) 

o 

entropy  production  (J  mol-1  K  ') 

er 

rate  of  entropy  production  (J  mol-1  K  1  s-1) 

electrolyte  conductivity  (S  m-1) 

Qd 

collision  integral  based  on  the  Lennard-Jones 

potential 

Z 

tortuosity 

Subscripts 

a 

anode 

A 

A  specie 

B 

B  specie 

c 

cathode 

cv 

control  volume 

h2 

hydrogen 

h2o 

water  vapour 

i 

mixture  component 

k 

Knudsen  diffusion 

n2 

nitrogen 

02 

oxygen 

p 

products 

R 

reactants 

(eft) 

effective 

Superscripts 

I 

inlet  condition 

Very  often,  the  Tafel  equation  or  the  linear  current- 
voltage  relation  is  used  to  replace  the  Butler- Volmer  equa¬ 
tion  in  the  calculation  of  activation  overpotential.  Or,  the 
Knudsen  diffusion  is  neglected  in  the  gas  diffusion  through 
the  porous  electrodes.  To  avoid  the  ambiguity  of  the  suit¬ 
ability  of  such  models  used  in  the  calculation  of  cell 
polarization,  the  Butler- Volmer  equation  is  retained  and 
both  the  ordinary  and  Knudsen  diffusions  are  considered 
in  the  proposed  polarization  model.  The  development  of  this 
model  forms  the  basis  of  the  objective  of  work  presented 
here.  In  addition,  an  analysis  is  conducted  to  show  the 
sensitivity  of  the  cell  voltage  drop  on  changes  in  the  thick¬ 
ness  of  respective  cell  components  under  fixed  exchange 
current  densities. 


2.  Thermodynamic  principles 

A  fuel  cell  is  just  like  any  other  energy  conversion  device/ 
system  in  that  its  performance  can  be  assessed  by  the 
principles  of  thermodynamics  [3,4].  Treating  the  fuel  cell 
as  a  control  volume  and  applying  the  steady  flow  energy 
equation  with  the  assumption  of  negligible  change  of  kinetic 
energy  and  potential  energy  gives 

Gcv  -  Wcv + y^ih,h,)R  -  y^(/iA)p = o 


(1) 
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The  entropy  balance  for  the  control  volume  can  be  written  as 

~  E(^')r  “  +  ^cv  =  0  (2) 

i  i 

Combining  Eqs.  (1)  and  (2)  and  assuming  that  the  electro¬ 
chemical  reaction  is  based  on  unit  molar  flow  rate  of  the  fuel 
(for  example  H2),  yields 


Wcv 


^2(ni8i)p  -  E^'^R 

i  i 


Tocv  =  0  (3) 


where  gt{T,pi )  =  ht(T)  -  TSi(T,Pi)  or  gi{T,Pi)  =  g°t  {T)+ 
RT  In  (p*//?o)  Since  there  is  no  mechanical  work  involved  in 
the  fuel  cell,  Wcv  is  naturally  the  electrical  work  done  by  the 
fuel  cell.  That  is 


We,  =  zFE 


E('^)p _  E  ("«'&)  r 


ToCY  (4) 


zFE  =  -A G°  - 


-  To cv 


where  — AG°  =  RT\nK(T). 

For  a  hydrogen  fuel  cell  with  known  pressures  of  reactants 
and  product 


where  [>  is  the  transfer  coefficient  and  i0  the  exchange 
current  density.  The  transfer  coefficient  is  considered  to 
be  the  fraction  of  the  change  in  polarization  that  leads  to  a 
change  in  the  reaction  rate  constant,  and  its  value  is  usually 
0.5  for  the  fuel  cell  application.  The  exchange  current 
density  is  the  forward  and  reverse  electrode  reaction  rate 
at  the  equilibrium  potential.  A  high  exchange  current  density 
means  that  a  high  electrochemical  reaction  rate  and  good 
fuel  cell  performance  can  be  expected. 

When  [>  =  0.5,  Eq.  (7)  can  be  expressed  as  follows: 


i  =  Hq  sinh 


(  n?Fh  acA 
V  2 RT  ) 


(8) 


Then, 


h  Act 


2  RT 
neF 


sinh  1 


(9) 


Two  special  cases  of  Eq.  (7)  are  available:  (1)  high  activation 
polarization;  (2)  low  activation  polarization. 

Under  high  activation  polarization,  the  first  exponential 
term  in  Eq.  (7)  will  be  much  greater  than  the  second  term, 
and  the  latter  can  be  excluded  from  the  equation.  Rearran¬ 
ging  the  simplified  equation  yields 


h  Act 


in  i'o  + 


(10) 


E 


RT 

—  in  K 
IF 


RT 
—  in 
4  F 


(Ph2o)Vo 

(Ph2)2P1o2. 


GcvT 

2  F 


(6) 


where  ocvT /IF  =  t)Aci  +  //()hm  +  //Conc. 

The  first  term  on  the  right-hand  side  of  Eq.  (6)  shows  the 
effect  of  the  temperature  on  the  cell  voltage  while  the  second 
term  shows  the  effect  of  the  pressures  of  the  reactants  and 
product  on  cell  voltage.  The  irreversibility  in  the  voltage 
drop  is  accounted  by  the  third  term  of  the  equation  and  is 
expressed  by  the  activation,  ohmic  and  concentration  over¬ 
potentials. 


2.1.  Activation  polarization 


Chemical  reactions,  including  electrochemical  reactions, 
involve  energy  banders  which  must  be  overcome  by  the 
reacting  species.  This  energy  bander  is  called  the  ‘activation 
energy’  and  results  in  activation  or  charge-transfer  polariza¬ 
tion,  which  is  due  to  the  transfer  of  charges  between  the 
electronic  and  the  ionic  conductors.  The  activation  polar¬ 
ization  may  be  regarded  as  the  extra  potential  necessary  to 
overcome  the  energy  barrier  of  the  rate-determining  step  of 
the  reaction  to  a  value  such  that  the  electrode  reaction 
proceeds  at  a  desired  rate. 

Activation  polarization  is  normally  expressed  by  the  well- 
known  Butler- Volmer  equation  [5] 


i  =  ?'o  jexp( 


»cT>iAcl\ 

RT  ) 


exp 


»eT>/Actl  1 

RT  \  J 


(7) 


which  is  the  well-known  Tafel  equation. 

Under  low  activation  polarization,  the  term  [l}]AciF/RT  in 
Eq.  (7)  will  be  much  less  than  unity  and  the  exponential 
terms  can  be  expanded  as  a  Taylor  series.  Neglecting  the 
terms  with  order  higher  than  unity  yields 


Act 


RT 

- i 

neFi0 


(11) 


which  is  the  well-known  linear  current-potential  relation. 

It  is  common  to  question  the  range  of  activation  polariza¬ 
tion  to  which  the  Tafel  equation  and  linear  current-potential 
relation  apply.  Thus,  a  plot  to  display  the  errors  of  the 
simplified  activation  models  relative  to  the  Butler- Volmer 
equation  will  be  useful.  The  errors  of  the  Tafel  equation  and 
the  lineal'  current-potential  relation  relative  to  the  Butler- 
Volmer  equation  at  1073  K  are  given  in  Fig.  1.  It  clearly 
shows  the  applicable  ranges  of  these  simplified  models.  For 
example,  if  5%  of  error  is  allowed,  which  corresponds  to 
Act  >  0.28  V,  Tafel  equation  can  be  used;  however,  if 
17 Act  <  0.1  V,  the  linear  relation  can  be  used.  For  a  SOFC, 
the  former  is  normally  applied  to  the  cathode  activation, 
while  the  latter  is  normally  applied  to  the  anode  activation. 


2.2.  Ohmic  polarization 


Ohmic  losses  occur  because  of  resistance  to  the  flow  of 
ions  in  the  electrolyte  and  resistance  to  flow  of  electrons 
through  the  electrode  materials.  The  dominant  ohmic  losses, 
through  the  electrolyte,  can  be  reduced  by  decreasing  the 
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Fig.  1.  Error  of  Tafel  and  linear  relation  results  relative  to  Butler- Volmer 
equation  at  1073  K. 

distance  of  electrode  separation  and  enhancing  the  ionic 
conductivity  of  the  electrolyte.  Because  the  ionic  flow  in  the 
electrolyte  obeys  Ohm’s  law,  the  ohmic  losses  can  be 
expressed  by  the  equation 

7  Ohm  =  iRe  (12) 

where  Re  is  the  ionic  resistance. 


(13) 


For  straight  and  round  pores,  the  diffusion  coefficient  of  the 
A  component  becomes 


(14) 


The  mean  pore  radius  can  be  evaluated  from  the  bulk 
density,  the  surface  area  of  the  porous  solid,  and  the  porosity 
by  the  expression 


2s 

r  = - 

SaPb 


(15) 


To  account  for  the  tortuous  path  of  the  molecule  rather  than 
along  the  radial  direction  and  the  porosity  of  the  solid 
electrode  for  the  fact  that  diffusion  occurs  only  in  the  pore 
space,  an  effective  Knudsen  diffusion  coefficient  may  be 
used,  i.e. 

^Ak(eff)  =  U\k  (16) 


Using  the  Chapman-Enskog  theory  of  prediction,  the  binary 
ordinary  diffusion  coefficient  in  the  gas  phase  can  be 
calculated  [7]  as  follows: 


2.3.  Concentration  polarization 

The  rate  of  mass  transport  to  the  reaction  sites  in  a  porous 
electrode  of  a  SOFC  can  be  described  by  the  diffusion  of 
gases  in  the  pores.  The  gases  have  to  diffuse  through  the  gas- 
filled  pores  of  the  electrode  in  order  to  reach  the  reaction 
sites.  Whenever  this  happens,  it  is  possible  for  the  concen¬ 
tration  polarization  to  be  significant,  since  when  the  current 
is  being  drawn  the  gas  partial  pressure  at  the  reaction  sites 
will  be  less  than  that  in  the  bulk  of  the  gas  stream.  Hence, 
decrease  of  gas  concentration  in  the  gas-filled  pores  of  the 
electrode  may  give  rise  to  severe  polarization  and  lead  to 
limiting  current.  Normal  concentration  polarization  does  not 
give  rise  to  excessive  voltage  losses  until  the  current  den¬ 
sities  approach  the  limiting  current. 

Diffusion  through  the  porous  material  is  typically 
described  by  either  ordinary  or  Knudsen  diffusion  and 
has  been  found  to  play  an  important  role  in  catalytic  reac¬ 
tion.  Ordinary  diffusion  occurs  when  the  pore  diameter  of 
the  material  is  large  in  comparison  to  the  mean  free  path  of 
the  gas  molecules.  Molecular  transport  through  pores  which 
are  small  in  comparison  to  the  mean  free  path  of  the  gas  is 
regarded  as  a  Knudsen  type  of  diffusion.  For  Knudsen 
diffusion,  molecules  collide  more  frequently  with  the  pore 
walls  than  with  other  molecules.  Upon  collision,  the  atoms 
are  instantly  adsorbed  on  to  the  surface  and  are  then 
desorbed  in  a  diffusive  manner.  As  a  result  of  frequent 
collisions  with  the  wall  of  the  pore,  the  transport  of  the 
molecule  is  impeded.  The  Knudsen  diffusion  coefficient  can 
be  predicted  using  kinetic  theory  by  relating  the  diameter  of 
the  pore  and  the  mean  free  path  of  the  gas  [6] 


/  1  1  \1/2  r3/2 

Dab  =  0.0018583  —  +  —  ,  - -  (17) 

\Ma  MbJ  pa2ABQDAB 

where  p  is  the  total  pressure  (in  atm),  erAB  =  (erA  +  oB)/2  is 
the  collision  diameter  (in  A)  and  £2DAB  is  the  collision 
integral  based  on  the  Lennard-Jones  potential  which  can 
be  obtained  from  eAB.  eAB  is  the  energy  of  molecular 
interaction  in  ergs  (eAB /k  =  \J (e A  / k)  (eB  / k) ,  K). 

Similar  to  the  Knudsen  diffusion,  the  effective  diffusion 
coefficient  for  ordinary  diffusion  can  be  described  by  the 
following  equation: 


^AB(eff) 


(18) 


Both  ordinary  diffusion  and  Knudsen  diffusion  may  occur 
simultaneously.  For  self-diffusion  or  equimolar  counter 
transfer,  the  effective  diffusion  coefficient  can  be  written  as 


A\(eff)  ^AB(eff)  A\k(eff) 

and 

^A(eff)  E  VU\B  DAk 

In  this  context,  the  combined  effect  of  both  types  of  diffusion 
on  concentration  polarization  is  considered.  If  it  is  assumed 
that  the  pores  of  different  sizes  are  highly  interlinked  and 
randomly  arranged,  then  it  can  be  assumed  that  the  path  of 
any  molecule  transfers  through  the  solid  electrode  should 
involve  essentially  the  same  flow  resistance  as  the  path  of 
any  other  molecules. 


(19b) 
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2.3.1.  Anode  —  equimolar  counter-current  diffusion  of 
reactant  and  product 

On  the  anode  side,  the  chemical  reaction  occurs  at  the  end 
of  the  diffusion  path.  The  flux  ratios  at  steady-state  are 
determined  by  the  stoichiometry  of  the  reaction.  Therefore, 

J  =  Jh2+  Jh2o  =  0  (20) 

Since  the  mass  of  hydrogen  diffusion  in  one  direction  does 
not  necessarily  equal  the  mass  diffusion  in  the  other  direc¬ 
tion,  a  definite  net  mass  flow  in  one  direction  is  present.  The 
equations  of  diffusion  are 

Jh2  =  ~ OH2(eff)VCn2  +  Ch2v  (21a) 

7h20  =  “ ^H20(eff)  VCh20  +  C\i2oV  (21b) 

where 

1  _  i  1 

Al2(eff)  Al2,k(eff)  Al2-H20(eff) 

1  _  1  1 

Ai20(eff)  Al20,k(eff)  Al2-H20(eff) 

For  steady  state  conditions,  VCh2  =  —  V C’h2o ,  Eqs.  (21a) 
and  (21b)  become 

Jh2  =  -Ai2(cff)VCH2  +  Ch2v  (22a) 

Jn2o  =  ^HoOfeff)  VCh2  +  Ch2()V  (22b) 

For  equimolar  counter-current  mass  transfer,  Jh2  =  —./mo 

—  (An(eff)  —  E)H2o(eff))VCH2  +  v(Ch2  +  Ch2o)  =  0  (23a) 


or 


Z^H2  (eff )  “  E>H2o(cff) 


VC, 


h2 


(23b) 


where  Ca  =  CH2  +  CH2o- 

Substituting  into  Eqs.  (21a)  and  (21b),  gives 


Jh2  =  ~ 

Jh2o  =  ~ 


Ch2 


C, 
Ch2o 

~c7 


Al2(eff)  +  I  Al20(eff) 


Ai2(eff)  + 


Ch2 

c7 


D 


H20(eff) 


VCh2  (24a) 
VCh2o 

(24b) 


Jh2  =  -Al(cIT)VCh2 

Jh20  =  ZZa(elT)  VCjj20 

where 


/-^a(eff) 


Ch2 


JD H2(eff)  +  J  //h20(c1T) 


or 


//a(eff) 


(~) 


DH2(eff)  +  (  Y7"  )  //[-liOfclT) 
P  a 


(25a) 

(25b) 


Since  the  difference  between  Z)H,(eff)  and  Z)mo(eff)  is  not 
very  large,  when  the  hydrogen  partial  pressure  does 
not  change  critically,  -Da(eff)  may  be  considered  as  a 
constant. 

For  one-dimensional  diffusion, 


Jr2 


■Z-^a(eff) 


dCm 

~dT 


(26) 


Because  /H,  =  —(i/2 F)  and  dCH2  =  (d pujRT),  so, 


dPn2 


RTi 

2FZ)a(eff) 


dl 


(27a) 


Integrating  this  equation  with  the  boundary  condition 
Pn2  =  Fh,  at  anode  surface  yields 


RTi 


Pn2  ~Pu2  = 


2  FD 


a(eff) 


(27b) 


Therefore,  in  terms  of  the  current  density,  i,  partial  pressure 
of  hydrogen,  pu2.  and  water  vapour  partial  pressure,  pn2o,  at 
reaction  sites 


RTL 


P  h2  =  Fh, 


2  FD 


a(eff) 


!  ,  RTI,  . 

P  h2o  —  Fh2o  +  9J7n  1 

2FL»a(eff) 

where 

i  £.  t  i  + 


//h2(c1T)  £a  \Dh2, k  At2-H2o 


1 


i  +  i 


//H20(eff)  £a  \Al20,k  Al2- H20 


(28a) 

(28b) 


2.3.2.  Cathode  —  self-diffusion 

On  the  cathode  side,  the  flux  of  nitrogen  is  zero.  Thus  [5], 

Jo2  =  —  DCV  Cq,  +  Xo2do2Jo2  (29) 

where 


<5o2 


//o2,k(eff) 

A>,,k(eff)  +  77o2-N2(eff) 


For  one-dimensional  diffusion, 


Jo2  —  — 77o2(eff)  ^ 2  +  Xo2502Jo2 

Because  Jq2  =  —i/4F  and  dCo2  =  d pQjRT,  so 
dpo2  RTi 


Pc  ~  <\hP(h  4FT>o2(eff)Pc 


dl 


(30) 


(31a) 


Integrating  this  equation  with  the  boundary  condition 
po2  =  p\h  at  cathode  surface  yields 


( (Pc/tch)  -Po2\  _  <fhRTi 
\(Pc/$02)  -Po2J  4FD02(eff)Pc  c 


(31b) 
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Therefore,  the  partial  pressure  of  oxygen  at  the  cathode 
reaction  sites  is 


P  o2  = 


Pc 

<5o2 


where 


P\h 


(  Sq2RTUc  \ 
\AFDc(c[t)pcJ 


1 

^02(eff) 


Sc 


and 


(32) 


A:(eff)  —  -Do2(eff) 


3.  Overall  cell  potential 


Polarization  is  a  very  important  parameter  in  the  analysis 
of  fuel  cell  performance.  For  calculating  the  relationship 
between  voltage  and  current  density,  the  following  general 
expression  can  be  used: 

E(i)  Eo  l?Ohm  ^Act.a  ^Act,c  Pc  one, a  ^7Conc,c 

(33a) 

where 


RT 

E0  =  —  In  K 
IF 


{(Afp'oJ 


P  Ohm  —  ‘Ec 


P  Act,  a 


2  RT 
neF 


sinh  1 


1  Act,c 


2  RT 
neF 


sinh  1 


RT  (l-(^r/2F)(Za/Da(eff)^2)i 
2  F  (l  +  (RT/2F)(la/D,(eff)plH20)i 


(33b) 

(33c) 

(33d) 

(33e) 

(33f) 


Table  1 


Calculation  based  parameters  for  solid  oxide  fuel  cell 


Parameter 

Value 

Fuel  cell  operating  temperature  ( T)  (K) 

1073 

Fuel  cell  operating  pressure  (p )  (atm) 

1 

Hydrogen  humidified  temperature  (TH)  (K) 

298 

Electrolyte  conductivity  (cre)  (S  m  ') 

10 

Electron  transferred  per  reacting  (nf) 

1 

Transfer  coefficient  (f) 

0.5 

Anode  exchange  current  density  (iQf)  (A  m~2) 

5300 

Cathode  exchange  current  density  (ioc)  (A  m~2) 

2000 

Cathode  thickness  (/c)  (m) 

50  x  10~6 

Porosity  (a)  (%) 

30 

Tortuosity  (c) 

6 

Average  pore  radius  (?)  (m) 

0.5  x  10~6 

For  electrolyte-supported  cell 

Electrolyte  thickness  (/e)  (m) 

500  x  10~6 

Anode  thickness  (If)  (m) 

50  x  10~6 

For  anode-supported  cell 

Electrolyte  thickness  (/e)  (m) 

40  x  10~6 

Anode  thickness  (If)  (m) 

750  x  10~6 

Current  density  (A/m2) 


Anode  —A— Cathode 


Fig.  2.  Calculated  activation  polarization  of  solid  oxide  fuel  cell. 


nCo 


RT 

= - In 

4  F 


(Pc/$o2)  -  {(Pc/$o2)  -  p1q2)  exp[(RT/4F)(So2lc/Dc(eff)Pc)i 


Po2 


(33g) 


Some  experimental  data  extracted  from  the  literature  [8-10] 
is  presented  in  Table  1.  Note  that  these  performance-related 
data  are  extracted  from  literature  which  focuses  on  the 
development  of  intermediate  temperature  SOFCs  that  oper¬ 
ate  at  a  temperature  around  800°C  [8-10].  Based  on  these 
data,  the  respective  polarization  curves  can  be  generated. 

3.1.  Activation  polarization 

Calculated  anode  and  cathode  activation  polarization 
curves  are  shown  in  Fig.  2.  The  cathode  polarization  is 
obviously  higher  than  that  of  the  anode  due  to  its  lower 


exchange  current  density.  The  activation  polarization 
increases  steeply  at  low  current  density  and  gradually  at 
high  current  density.  The  convex-up  curvature  of  polariza¬ 
tion  is  the  cause  for  the  concave  voltage  versus  current 
density  trace.  It  should  be  noted  that  the  activation  loss  in  the 
present  study  was  assumed  to  occur  at  the  electrode/elec¬ 
trolyte  interface  (or  two-dimensional).  In  an  actual  fuel  cell, 
the  electrochemical  reaction  occurs  at  electrochemical 
active  sites  (i.e.  three-dimensional)  where  the  reactive 
gas,  the  electronic  conductor  and  the  ionic  conductor  are 
present  at  the  same  time.  In  other  words,  the  thickness  of  the 
electrode  does  have  an  effect  of  the  activation  loss,  which  is 
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Current  density  (A/m2) 


—♦—Anode  —A— Cathode 


Fig.  3.  Calculated  concentration  polarization  of  anode-supported  solid 
oxide  fuel  cell. 


Fig.  4.  Calculated  concentration  polarization  of  electrolyte-supported 
solid  oxide  fuel  cell. 


ignored  in  the  present  study.  The  present  model,  however, 
can  be  extended  to  cover  the  micro  level  of  electrode  kinetics 
based  on  the  work  of  Costamagna  et  al.  [11]. 

3.2.  Concentration  polarization 

Anode  concentration  polarization  and  cathode  concentra¬ 
tion  polarization  are  calculated  using  Eqs.  (33f)  and  (33g), 
respectively.  The  calculated  concentration  polarization  of  an 
anode-supported  SOFC  is  shown  in  Fig.  3.  The  anode 
concentration  polarization  is  much  higher  than  that  of  the 
cathode  because  of  the  thick  anode  used.  It  should  be 
noted  that  the  polarization  curve  is  convex  at  low  current 
density.  At  high  current  density,  the  polarization  curve 
exhibits  a  concave  curvature.  Whether  or  not  a  concave 
curvature  is  observed  at  high  current  density  (in  the  cell 
voltage-current  density  trace)  depends  really  on  the  magni¬ 
tude  of  the  ohmic  and  activation  polarization.  For  the  data 
in  Table  1,  the  cell  voltage  drops  to  zero  when  the  current 
density  reaches  26,600  A  m-2.  Therefore,  this  phenomenon 
(concave  curvature)  cannot  be  seen.  Actually,  the  selected 
parameters  of  exchange  current  densities,  electrolyte  con¬ 
ductivity  and  thickness  are  amongst  the  highest  values 
obtained  from  the  literature.  It  can  be  concluded  that  the 
anode  limiting  current  density  is  not  reached  under  normal 
conditions. 

The  calculated  concentration  polarization  of  an  electro¬ 
lyte-supported  SOFC  is  given  in  Fig.  4.  The  curve  is  slightly 
convex  over  the  whole  range  of  current  density.  Because  its 
magnitude  is  small,  it  will  not  influence  the  shape  of  the  cell 
voltage  versus  current  density  trace. 

3.3.  Cell  voltage  and  power  density 

The  cell  voltage  and  power  densities  are  calculated  for  the 
whole  range  of  fuel  cell  operations  from  i  =  0  to  V  =  0. 

The  calculated  cell  voltage,  polarization  and  power  den¬ 
sity  versus  current  density  trace  of  the  electrolyte- supported 


cell  are  presented  in  Fig.  5.  The  power  density  reaches  a 
peak  value  when  the  current  density  is  5500  A  in  2.  At  this 
point,  the  voltage  is  0.519  V  and  the  power  density  is 
2834  W  m  2.  Normally,  the  fuel  cells  are  designed  to  oper¬ 
ate  slightly  left  to  this  point  for  a  good  compromise  between 
cell  efficiency  and  low  capital  cost  and  the  consideration  of 
operational  stability.  Under  this  condition,  the  maximum 
loss  is  due  to  the  electrolyte  ohmic  polarization,  and  the 
cathode  and  anode  activation  losses  are  also  considerable. 
The  concentration  losses  are  relatively  small  in  this  case. 
When  the  current  density  reaches  1 1 ,500  A  m  2,  the  cell 
voltage  drops  to  zero.  From  this  study,  one  may  identify  the 
right  polarization  to  be  focused  on  for  improvement.  The 
electrode- supported  cells,  which  have  thinner  electrolyte, 


— ♦—  Ctel I  vol  t  age 
— B—  Chni  c  pol  ar  i  sat  i  on 
— a —  Axxte  act  i  vat  i  on  pol  ar  i  sat  i  on 
-*-Cht  hode  act  i  vati  on  pol  ar i  sati  on 
— finode  concent  rat  i  on  pol  ari  sat  i  on 
— e—  Cht  hode  concent  r  at  i  on  pol  ar  i  sat  i  on 
■  Ftwer  densi  t  y 

Fig.  5.  Calculated  cell  voltage,  polarization  and  power  density  vs.  current 
density  trace  of  electrolyte-supported  cells. 
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Cell  voltage  —a— Ohmic  polarisation 

—a—  Anode  activation  polarisation  — x—  Cathode  activation  polarisation 

— Anode  concentration  polarisation  — e—  Cathode  concentration  polarisation 
— Power  density 


Fig.  6.  Calculated  cell  voltage,  polarization  and  power  density  vs.  current 
density  trace  of  anode-supported  cells. 


could  have  better  performance  than  that  of  electrolyte- 
supported  cells. 

The  calculated  cell  voltage,  and  power  density  versus 
current  density  trace  of  an  anode-supported  cell  are  given  in 
Fig.  6.  As  discussed  above,  it  can  be  concluded  that  the 
concave  curvature  of  the  cell  voltage  is  due  to  both  activation 
and  concentration  polarization.  The  power  density  reaches  a 
peak  value  when  the  current  density  is  12,000  A  m  2.  At  this 
point,  the  voltage  is  0.41  V  and  the  power  density  is 
4900  W  m  2.  Under  such  condition,  the  maximum  loss  is 
due  to  the  cathode  activation  polarization  and  the  anode 
activation  and  concentration  losses  are  also  significant.  The 
electrolyte  ohmic  and  cathode  concentration  losses  are 
relatively  small.  When  the  current  density  reaches 
26,600  A  m  2,  the  voltage  drops  to  zero.  Therefore,  the 
electrode  material-related  factors  such  as  the  catalytic  acti¬ 
vity  and  reaction  sites,  which  determine  the  magnitude 
of  the  exchange  current  density,  should  be  targetted  for 
improvement. 

Comparing  the  structure  of  anode- supported  and  electro¬ 
lyte-supported  cells,  the  only  difference  is  the  thickness 
of  the  anode  and  the  electrolyte.  Therefore,  in  general, 


It  should  be  noted,  at  this  point,  that  the  effects  of  the 
porosity,  tortousity  and  pore  size  of  the  electrodes  on  the  fuel 
cell  performance  are  not  considered  because  the  present 
model  does  not  examine  the  effects  of  these  parameters  on 
the  effective  electrochemical  reaction  sites.  Hence,  it  is  not 
appropriate  to  show  their  effects  on  the  concentration  over¬ 
potential  without  considering  their  impacts  on  the  electro¬ 
chemical  reaction.  As  mentioned  in  Section  3.1,  once  the 
microelectrode  kinetic  model  is  included,  the  interaction 
between  the  activation  and  concentration  overpotentials  due 
to  the  change  of  porosity,  tortousity  and  pore  size  will  be 
meaningful. 


4.  Sensitivity  analysis 


The  sensitivity  study  of  the  proposed  model  is  performed 
by  defining  the  dimensionless  sensitivity  of  the  cell  voltage 
subject  to  the  change  of  any  parameter  of  interest  (P,) 


_  dE/E 
~  dPi/Pi 


(34) 


For  a  change  in  anode  thickness 
dE/dla 


SL  = 


E/l  a 


(35a) 


where 


dE 

dla 


+ 


RT 


X 


2  E  )  \T>a(eff) 


Pu2  o  (RT/2E)  (i/T)a(e|-fj)/a 


Ph2o  +  (f?T  /  2F)  (i/ Da(et-f  j )  /a 


(35b) 


For  a  change  in  electrolyte  thickness, 

c  _  dE/ dle 
k~  E/lc 

where 
dE  _  i 
dle  OJI 

For  a  change  in  cathode  thickness, 

c  _  dE/  dlc 
k~  E/lc 

where 


(36a) 


(36b) 


(37a) 


dE  _  XRT\2(  1  -  {PoJPc)^o2 )  exp[(/?r/4Fpc)(fi02/c/T>c(eff))/](//Dc(eff)) 

die  V4iV  (Pc/<5 02 )  -  ((Pc/<5 02 )  -  Po2)  ^M(ET/4FPe)(So2lc/De(eff))i\ 


it  can  be  concluded  that  the  performance  of  anode-supported 
cells  should  be  better  than  that  of  electrolyte- supported 
cells. 


The  sensitivity  of  the  cell  voltage  due  to  the  change  of 
respective  cell  component  thickness  at  different  current 
densities  is  shown  in  Figs.  7-11.  When  all  three  cell 
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(c)  I  (A/m2) 


Fig.  7.  (a)  Sensitivity  of  cell  voltage  due  to  anode  thickness  (anode/ 
cathode/electrolyte:  500  (im).  (b)  Sensitivity  of  cell  voltage  due  to 
electrolyte  thickness  (anode/cathode/electrolyte:  500  |tm).  (c)  Sensitivity 
of  cell  voltage  due  to  cathode  thickness  (anode/cathode/electrolyte: 
500  |im). 


Fig.  8.  (a)  Sensitivity  of  cell  voltage  due  to  anode  thickness  (anode/ 
cathode/electrolyte:  250  (im).  (b)  Sensitivity  of  cell  voltage  due  to 
electrolyte  thickness  (anode/cathode/electrolyte:  250  |tm).  (c)  Sensitivity 
of  cell  voltage  due  to  cathode  thickness  (anode/cathode/electrolyte: 
250  |im). 
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Fig.  9.  (a)  Sensitivity  of  cell  voltage  due  to  anode  thickness  (anode: 
750  (.im;  cathode:  50  (.tm;  electrolyte:  40  |im).  (b)  Sensitivity  of  cell 
voltage  due  to  electrolyte  thickness  (anode:  750  (tm;  cathode:  50  |im; 
electrolyte:  40  (.im).  (c)  Sensitivity  of  cell  voltage  due  to  cathode  thickness 
(anode:  750  |im;  cathode:  50  (im;  electrolyte:  40  (im). 


components,  namely,  anode,  cathode  and  electrolyte,  have 
the  same  thickness  of  500  pm  (Fig.  7a-c),  S/e  is  the  largest, 
then  S[c  and  finally  S/a .  Decreasing  all  three  thicknesses  to 
250  pm  (Fig.  8a-c),  in  general,  results  in  reduced  sensitiv¬ 
ities  and  an  improved  operating  current  density  range. 
Intuitively,  these  sensitivity  profiles  show  that  the  use  of 
electrolyte  as  the  fuel  cell  support  is  not  as  good  as  that  using 
either  anode  or  cathode  as  the  support.  The  sensitivity 
profiles  of  anode-supported  and  cathode-supported  cells 
are  compared  in  Figs.  9  and  10,  respectively,  The  results 
show  that  the  anode- supported  cell  is  superior  to  the  cath¬ 
ode-supported  cell  in  terms  of  operating  current  density 
range  and,  hence,  higher  power  density  can  be  achieved. 
It  is  obvious  that  anode  concentration  overpotential  is  the 
cause  of  limiting  current  for  the  former,  whereas  cathode 
concentration  oveipotential  is  the  cause  of  limiting  current 
for  the  latter.  To  compensate  for  the  serious  cathode  con¬ 
centration  oveipotential  loss  in  the  cathode- supported  cell, 
the  cathode  compartment  pressure  is  increased  from  atmos¬ 
pheric  pressure  to  4  bar.  In  general,  all  three  sensitivities  are 
reduced  with  improved  operating  current  density  range,  see 
Fig.  1  la-c.  Compared  with  an  anode- supported  cell  operat¬ 
ing  at  atmospheric  condition,  however,  the  cathode-sup- 
ported  cell  is  no  way  better  in  terms  of  performance. 
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Fig.  10.  (a)  Sensitivity  of  cell  voltage  due  to  anode  thickness  (anode: 
50  pm;  cathode:  750  pm;  electrolyte:  40  pm).  (b)  Sensitivity  of  cell 
voltage  due  to  electrolyte  thickness  (anode:  50  pm;  cathode:  750  |im; 
electrolyte:  40  pm).  (c)  Sensitivity  of  cell  voltage  due  to  cathode  thickness 
(anode:  50  pm;  cathode:  750  |tm;  electrolyte:  40  pm). 


5.  Summary  and  conclusions 

It  was  felt  that  a  generalized  cell  voltage-current  density 
model  could  be  important  and  useful  for  researchers  who  are 
working  on  SOFCs.  Much  effort  in  this  study  has  been 
devoted  to  the  development  of  this  generalized  model  based 
upon  thermodynamic  principles  whereby  the  irreversibility 
term  in  the  steady  flow  thermodynamic  system  is  manipu¬ 
lated  electrochemically  and  replaced  by  the  overall 
polarization. 

This  paper  points  out  the  applicable  range  of  the  Tafel 
equation  and  the  simplified  linear  current-potential  expres¬ 
sion  of  the  Butler- Volmer  model  which  are  easily  mis-used. 
If  5%  of  error  relative  to  the  Butler- Volmer  model  is  allowed 
for  both  cases,  the  activation  overpotential  for  the  Tafel 
equation  should  be  more  than  0.28  V  and  for  the  linear 
expression,  <0.1  V.  To  generalize  the  cell  voltage-current 
density  model,  the  original  Butler- Volmer  expression  is 
used  in  the  calculation  of  activation  polarization  in  conjunc¬ 
tion  with  the  consideration  of  both  ordinary  and  Knudsen 
diffusion  in  the  calculation  of  concentration  polarization. 

To  develop  a  high-performance  SOFC,  the  selection  of  a 
particular  cell  component  as  the  support/substrate  of  the  fuel 
cell  can  be  the  determining  factor.  Our  study  shows  that  the 


Fig.  11.  (a)  Sensitivity  of  cell  voltage  due  to  anode  thickness  (anode: 
50  pm;  cathode:  750  pm;  electrolyte:  40  pm;  pc :  4  bar),  (b)  Sensitivity  of 
cell  voltage  due  to  electrolyte  thickness  (anode:  50  pm;  cathode:  750  pm; 
electrolyte:  40  pin;  pc:  4  bar),  (c)  Sensitivity  of  cell  voltage  due  to  cathode 
thickness  (anode:  50  pm;  cathode:  750  pm;  electrolyte:  40  pm;  pc :  4  bar). 


performance  of  an  anode-supported  cell  is  superior  to  that  of 
an  electrolyte-supported  cell  or  a  cathode-supported  cell  for 
the  same  materials  used,  same  electrode  kinetics  and  same 
operational  constraints.  Results  show  that  when  all  three 
components  are  of  the  same  thickness,  the  sensitivity  of  cell 
voltage  due  to  the  change  of  electrolyte  thickness  is  the 
highest,  then  sensitivity  of  cell  voltage  due  to  the  change  of 
cathode  thickness  and  finally  the  sensitivity  of  cell  voltage 
due  to  the  change  of  anode  thickness.  Reducing  all  three 
component  thicknesses  causes,  in  general,  reduced  sensiti¬ 
vity  strength  with  an  improved  operating  current  density 
range.  The  study  also  reveals  that  a  cathode-supported  cell 
is  no  way  better  than  an  anode- supported  cell  for  high 
performance  fuel  cell  design  even  under  elevated  pressure 
conditions  at  the  cathode  compartment  to  compensate  for 
the  serious  cathode  concentration  overpotential  loss. 
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